Feedback control on geometric phase in dissipative two-level systems 
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The effect of feedback on a two-level dissipative system is studied in this paper. The results show 
that it is possible to control the phase in the open system even if its state can not be manipulated 
from an arbitrary initial one to an arbitrary final one. The dependence of the geometric phase on 
the control parameters is calculated and discussed. 
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Control on a quantum system is a major task required 
for quantum information processing. Several approaches 
to reach this goal have been proposed in the past decade, 
which can be divided into open loop and closed loop 
problems, according to how the controls enter the dy- 
namics. In the open loop problems, the controls enter 
the dynamics through the system Hamiltonian. It affects 
the time evolution of the system state, but not its spec- 
trum, i.e., the eigenvalues of the target density matrix 
pf remain unchanged in the dynamics due to the uni- 
tarity of the evolution. In the closed loop problems, a 
feedback is required from the controller to enter the sys- 
tem based on measurement performed for the open sys- 
tem. This feedback control strategy can be traced back to 
1980's[ll[l,Q when it is used to explain the observation^ 
of subshot-noise fluctuations in an in-loop photon cur- 
rent. In these works, the authors did this using quan- 
tum Langevin equations and semiclassical techniques, the 
latter approach was made fully quantum-mechanical by 
Plimak[5j. For linear systems, all of these approaches and 
the trajectory approach 0,0,9 are equally easy to use to 
find analytical solutions. Nevertheless, the trajectory ap- 
proach has advantages that (1) it is applicable for quan- 
tum system with non-linear dynamics, and (2) it is very 
easy to consider the limit of Markovian (instantaneous) 
feedback, i.e., a master equation describing the uncon- 
ditional system dynamics with feedback may be possible 
to derive [y, 0, The Markovian feedback can be used 
to modify the stationary state in a two-level dissipative 
quantum system, and this feedback manifests interest- 
ing effects on the time-optimal control in the quantum 
system governed by the Lindbald master equation 

Geometric phases in quantum theory attracted great 
interest since Berry[l3] showed that the state of a quan- 
tum system acquires a purely geometric feature in ad- 
dition to the usual dynamical phase when it is varied 
slowly and eventually brought back to its initial form. 
The Berry phase has been extensively studied IllL 
and generalized in various directions [l3|, Il4l. llSl lla . Il7l |. 
such as geometric phases for mixed states for open 
systems [ig], and with a quantized field driving 13] ■ In 
a recent paper [Til, Sjoqvist calculated the geometric 
phase for a pair of entangled spins in a time-independent 
uniform magnetic field. This is an interesting develop- 
ment in holonomic quantum computation and it shows 



how the prior entanglements modify the Berry phase. 
This study was generalized [l^ to the case of spin pairs 
in a rotating magnetic field, which showed that the ge- 
ometric phase of the whole entangled bipartite system 
can be decomposed into a sum of geometric phases of 
the two subsystems, provided the evolution is cyclic. 
A renewed interest in geometric phenomena in quan- 
tum physics has been recently motivated by the pro- 
posal of using geometric phases for quantum comput- 
ing. Geometric phases depend only on the geometry of 
the path executed, and are therefore resilient to certain 
type of errors. The idea is to explore this inherent ro- 
bustness provided by the topological properties of some 
quantum systems as a means of constructing built-in 
fault tolerant quantum logic gates. Various strategies 
have been proposed to reach this goal, some of them 
making use of purely geometric evolution [20l [2ll . [2^ . 
Others make use of hybrid strategies that combine to- 
gether geometric and dynamical evolution (23l [2^. Sev- 
eral proposals for geometric quantum computations have 
been suggested and realized in different context^ includ- 
ing NMR experientsdll, ion traps [gE [2I [13, IH, HI], 
cavity QED^SQl, atomic ensembles |3l|, [34|, Josephson 
j miction |33|. anyonic system[34:] and quantum dotfssl. 

For open systems governed by the Lindblad master 
equation, it was shown that the controls can not fully 
compensate the effect of decoherence, indicating that the 
state of open systems can not be manipulated from an 
arbitrary initial state to an arbitrary final state. This 
gives rise to a question that can the geometric phase of 
such an open system be controlled? 

In this paper, we shall study the effect of feedback 
on the geometric phase in a dissipative two-level system 
governed by the Lindblad master equation. Consider an 
atom with two relevant levels {\g), |e)} and lower opera- 
tor (T_ = |(7)(e|. Let the atomic decay rate be 7 and let it 
be driven by a classical magnetic field B{t). Within the 
Markovian approximation for the system-environment 
couplings, the time evolution of the two- level system is 
described by the Lindblad master equation, 
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C{p) = i-f{Fa-pa+F^ - ^pcr^cr' 



Here a closed-loop control F is introduced, which is trig- 
gered immediately only after a detection click, namely a 
quantum jump occurs. This scheme was used to generate 
and protect entangled steady state in cavity QED system 
and the jump feedback Fa- can be understood as 

follows. The unitary operator F is applied only immedi- 
ately after a detection event, which is described by term 
a-pa+. Intuitively the stationary states depend on the 
feedback operator F. So, once the measurement prescrip- 
tion has been chosen, the freedom to design a feedback to 
produce a stationary state lies in the different choices for 
the feedback operator F. Although an enormous range 
of possibilities for F is allowed, even considering the lim- 
itations imposed by experimental constraints, we here 
choose (with the constraint FF^^ ~ 1) 

F = e'^-^ ^cosA+ '^ ' ^ sin A, (2) 

A 

where we denote A= I A I . In fact the feedback F written 



in this form covers all allowed possibilities. Writing the 
reduced density matrix 



we show after a simple algebra that 
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and 



C{p) = ^■{—^l^--^ia --{- + -■ a) 

sin 2 A, — > „ A, sin . 
• [cos2A-^ + — - — -ctxA+2^^ A]}. 



(5) 



Here it = (0,0,1), and the master equation can be 
rewritten as. 
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For an open system, its state in general is not pure and 
the evolution of the system is not unitary. For non- 
unitary evolution, the geometric phase can be calculated 
as follows. First, solve the eigenvalue problem for the 
reduced density matrix p[t) and obtain its eigenvalues 
Ek{t) as well as the corresponding eigenvectors \Ek{t))] 
Second, substitute Ek{t) and \Ek{t)) into [37], 



7g(r) = arg^ 



{Ek{t = 0)|£;fc(r))e"r^^'=<*^l^'=^*"'^* 



(7) 

where 7g is the geometric phase for the system undergo- 
ing nonunitary evolution[37|. In the following we shall 
choose 



|*(0)) =cos-|e)+sin-|5) 



(8) 



as the initial state for the dissipative two-level system, in 
terms of p = {px,Py,Pz) the initial state can be expressed 
as 

= sin 6*, = 0, = cos 6*. (9) 
The final state of the dissipative two-level system is then 
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with a and (j) defined by. 



pI+pI+ pI 



tan< 
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(11) 



The eigenvalues and corresponding eigenstates of the re- 
duced density matrix follows 



1 1 
- ± - 
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^pI+pI+pI (12) 



and 



\E+) = cos|e'^|e) + sin||g). 



\E\ 



sm ^e*"^|e) - cos 



2 ' ' 2 
respectively. It is easy to check that, 
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y^p2(o)+p2(o)+p2(o)^0, (14) 
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FIG. 1: (color online) Contour plots for the geometric phase 
as a function of A (vertical axis) and l3 (horizontal axis), u 
was chosen to be oj = 0.005/X-Bo, and [0,r] represents time 
interval of the system evolution, r = 27r/tj. In these plots, we 
set IJ.B0 = 1, and both A and (3 are in units of vr. (a) 7 — 
0.001, (b) 7 = 0.005, (c)7 = 0.01, (d) 7 = 0.05, (0)7 = 0.1, 
(f)7 = 0.5, (g)7 = 1 and (h)7 = 3. 



so the geometric phase jg reduces to 

{E+{t = 0)|i;+(t))e"/o^<^+^*^'^+^*^^''* 



7k = arg 



straightforward calculations show that 



(£;+(i)|||£;+W>=»cos2(f)|^, 



and 
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We can use these equations to perform numerical 
simulations of the geometric phase for the open sys- 
tem. In the numerical simulation, we choose B{t) = 
Bf) (cos 8 cos $, cos 6 sin $, sin Q) as the varying mag- 
netic field with <i> = uit. Without atomic decay, i.e.. 





FIG. 2: (color online) The same as Fig[T] but with the vertical 
axis is 7. (a)^ = 7r/4, (b)A = 7r/2, (c)yl = 37r/4, (d)A = tt. 
The other parameters chosen are the same as in Fig[T] 



7 = 0, the two-level system evolves freely, and its dy- 
namics is governed by idp/dt = [Hq,p\, where i/o is 
given in Eq. ([T|) . The instantaneous eigenstates of are 
1+) =cosfe**|e)-fsinf |-) = sin f e'*|e)-cos f 
with the corresponding eigenvalues e± = ±fiBo, respec- 
tively. For this system to evolve adiabatically, the adia- 
batic condition requires lu <ti (J^Bq. In the numerical sim- 
ulation, uj was chosen to be a; = O.OOS/iSo, and r — 27r/a;. 
To be specific, we set = Asmf3,Ay = Acosf3, and 
Az = 0, where A is a constant. So the feedback is char- 
acterized by (3 and A. The plots presented in Fig[l]are for 
the geometric phases acquired by the dissipative two-level 
system as a function of A and f3. For very small atomic 
decay rate 7^0 (FiglT]-(a)), the geometric phase ap- 
proaches a constant 7^ ~ (1 — cos 9) (in units of tt). As 7 
increases, the range of the geometric phase acquired by 
the dissipative system increases, implying that the geo- 
metric phase can be controlled even with large atomic 
decay rate 7. This is different from the control on quan- 
tum states, where the control can not fully compensate 
the decoherence. Figure [1] also shows that the geomet- 
ric phase is a periodic function of A, this can be under- 
stood by examining Eq.([21), where the feedback control 
is given. In addition to the above observation, we can 
find from figure [1] that the geometric phase is regular for 
small and large 7, while it is irregular for intermediate 
values of 7. The physics behind this feature is the fol- 
lowing. For very small 7, C{p) is negligible, hence Hq 
dominates over C{p) in the dynamics and the geomet- 
ric phase is mainly determined by Hq. When 7 is large 
enough such that C{p) dominates the dynamics, the geo- 
metric phase then comes from the dynamics governed by 
C{p). With a specific A — 7r/4, the geometric phase as 
a function of (3 and 7 is plotted in figure O In all these 
plots, we set Q = 9, thus the initial state is an eigenstate 
of the free Hamiltonian Hq. For A = 7r/2, the feedback is 
F ~ i{sinf}ax + cos fiay) whereas for A = tt, F becomes 
1, i.e., there is no feedback operating on the system. This 
can be found in Figl21-(d), where the geometric phase ac- 
quired is independent of f3. From figure [21-(a) and[2]-(c) 
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we find that figure [2]-(c) is exactly the same as figure 
[2]-(a) by replacing /3 by /3 + tt, indicating that the geo- 
metric phase remains unchanged with A ^ n ~ A and 
P ^ P + TT. This feature can be understood as follows. 
Recall that A = {Asm (3, Acos (3,0), F can be written as 
F — cos A + i{ax sin (3 + ay cos (3) sin A, it is clear that F 
remains unchanged by replacing A and (3 with t: — A and 
TT + /3, respectively. 

To sum up, in this paper, we have studied the effect 
of feedback on the geometric phase of a dissipative two- 
level system. The dependence of the phase on the feed- 



back parameters are calculated and discussed. The re- 
sults suggested that we can manipulated the phase by a 
properly designed feedback control. For small and large 
atomic dissipative rates with respect to the amplitude of 
the driving magnetic field ^Bq, the geometric phase is a 
periodic function of the feedback parameters, the physics 
behind these features is also presented. 
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